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Abstract 

Properties and temperature behavior of it and a bound states are studied in the 
framework of the nonlocal model with a separable interaction kernel based on the 



quark Dyson- Schwinger and the meson Bethe-Salpeter equations. M^(T), fw(T), 



M a (T) and r cr ^ 7r7r (T) are considered above and below the deconfinement and chiral 
restoration transitions. 



1 Introduction 

Understanding the behavior of matter under extreme conditions is nowadays a challenge in the 
physics of strong interactions. Different regions of the QCD phase diagram are an object of 
interest, and major theoretical and experimental efforts have been dedicated to the physics of 
relativistic heavy-ion collisions looking for signatures of the quark gluon plasma [QGP] [H 121 El- 
Restoration of symmetries and deconfinement are expected to occur at high-density and/or 
temperature. In this regard, the study of observables of pseudoscalar and scalar mesons is 



particularly important. Since the origin of these mesons is associated with the phenomena 
of spontaneous and explicit chiral symmetry breaking, its temperature behavior is expected 



Oh" 

to carry relevant signs of a possible restoration of symmetries. Usually, the restoration of 
chiral symmetry at high temperature is connected with the transition of hadron matter into 
quark-gluon plasma. 

Effective quark models are useful tools to explore the behavior of matter at temperatures. 
Nambu-Jona-Lasinio [NJL] [3] type models have been extensively used over the past years to 
describe low-energy features of hadrons and also to investigate restoration of chiral symmetry 
with temperature [Sj-jH]- 

This paper is devoted to investigation of the phase transition in hot matter and the tem- 
perature behavior of pseudoscalar as well as scalar mesons in the framework of the effective 
nonlocal model. This work is the continuation of . In [§] , a special separable form of effective 
gluon propagator is used in construction of the quark Dyson - Schwinger equation (DSE) and 
the Bethe - Salpeter equation (BSE) for bound states. Only pseudoscalar and vector mesons 
are considered in that paper. Here we concentrate on the properties of the scalar a meson at 
finite temperature. 

1 Talk at the Round Table Discussion "Searching for the mixed phase of strongly interacting matter 
at the JINR Nuclotron", JINR, Dubna, July 7 - 9, 2005. 
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2 Dyson - Schwinger equation with separable interaction 



The dressed quark propagator S(p) and meson Bethe - Salpeter (BS) amplitude T(p, P) are 
solutions of the DSE [E]-[TIj 

Sip)" 1 = ^p + m + ^J -0^9 2 Df(p - q)l,S{q) lv (1) 

and the BSE equation 

r(p, P) = ~ J -0^9 2 Dl% - q)^S(q + )T(q, P)%_) 7 „ (2) 

rr 

where 2 D^;(p — q) is an "effective gluon propagator", m is the current quark mass, P is the 
total momentum, and q± = q ± P/2. The form of equations (0) and (J2J) corresponds to the 
rainbow - ladder truncations of DSE and BSE. 

The simplest separable Ansatz g 2 D^(p — q) — > 5^ u D(p 2 , q 2 ,p- q) in a Feynman - like gauge 
is employed 

IV, g 2 ,p • g) = D f (p 2 )f (q 2 ) + D 1 f 1 (p 2 )(p ■ q)h{q 2 ) (3) 

This is a rank-2 interaction with two strength parameters D , D\ and the corresponding form 
factors fi{p 2 )- The choice for these quantities is constrained to the solution of the DSE for the 
quark propagator Q 

Sip)' 1 = ipA{p 2 ) + B(p 2 ) (4) 

with B(p 2 ) = mo + bfo(p 2 ) and A(p 2 ) = 1 + af\(p 2 ), where a and b are some constants. If there 
are no poles in the quark propagator S(p) for real timelike p 2 , then there is no physical quark 
mass shell and a meson cannot decay into quarks. The propagator is confining 3 if m 2 (p 2 ) ^ —p 2 
for real p 2 where the quark mass function is m(p 2 ) = B(p 2 )/A(p 2 ). It is found that the simple 
choice fi(p 2 ) = exp (—p 2 /A 2 ) leads to a reasonable description of tt and a meson properties. 
At the same time, the produced quark propagator is found to be confining and the infrared 
strength and shape of the quark amplitudes A(p 2 ) and B(p 2 ) are in qualitative agreement with 
the results of typical DSE studies. We use the exponential form factors as a minimal way to 
preserve these properties while realizing that the ultraviolet suppression is much greater than 
the power law fall-off (with logarithmic corrections) known from asymptotic QCD. The total 
number of model parameters is five. 

The extension of the separable model studies to the finite temperature case, T ^ 0, is 
systematically accomplished by a transcription of the Euclidean quark 4 - momentum via q — > 
q n = (uj n ,q), where u n = {2n + 1)itT are the discrete Matsubara frequencies. The effective qq 
interaction will automatically decrease with increasing T without the introduction of an explicit 
T-dependence which would require new parameters. 

The result of the DSE solution for the dressed quark propagator now becomes 

S-\p n ,T) = i^-pA{ pl T) + i l4 ou n C(pl T) + B(p 2 n , T) , (5) 

2 We use the Euclidean metric. 

3 A similar confining propagator without poles is used in a nonlocal model of NJL tvne |l()l [TT] . 
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where p 2 n = u 2 + p 2 . The solutions have the form B = mo + &(T)/ (p 2 ), A = 1 + a(T)/i(p^), 
and C — 1 + c(T)/i(p^), and the DSE becomes a set of three non-linear equations for 6(T), 
a(T), and c(T). The explicit form is 

<T) = ^T^J^f^p'll + a^f^d-^plT), (6) 
c(T) = '^/^^[l + ^KU^, (7) 

= ^r^/^/o^K + K^/o^)]^ 1 ^^), (8) 

n 

where d(p 2 ,T) is given by 

% 2 , 20 = p 2 A 2 (p 2 , T) + W 2 C 2 (p 2 , T) + P 2 (p 2 , T). (9) 

3 Pseudoscalar and scalar correlations 

With the separable interaction the allowed form of the solution of (J2J) for the pion BS amplitude 
fr 7r (g;P) is 

I\r(g; P) = 7s (^(P 2 ) + PF^P 2 )) f (q 2 ) . (10) 

The q dependence is described only by the first form factor fo{q 2 )- The second term fx of the 
interaction can contribute to BS amplitude only indirectly via the quark propagators. The pion 
BSE, (J2J), becomes a 2 x 2 matrix eigenvalue problem JC(P 2 )f = X(P 2 )f where the eigenvector 
is / = (E n , F n ). The kernel is 

MP 2 ) = "^tr. j tJ^/oV) [US(q + )t 3 S(q^)] , (11) 

where the it covariants are t = (ijs, 75 P) with t = (275, —75 P/2P 2 ). We note that the separable 
model produces the same q 2 shape for both amplitudes F n and E.,,; the shape is that of the 
quark amplitude B(q 2 ). Goldstone's theorem is preserved by the present separable model; in 
the chiral limit, whenever a nontrivial DSE solution for B(p 2 ) exists, there will be a massless 
71 solution to (fTTj) . 

The normalization condition for the pion BS amplitude can be expressed as 



2P » = 2Nc dF; J ^M^(?;-^)%+)r.(g;in%-)) 



(12) 

p2=K 2 =-Ml 



Here T(q;K) is the charge conjugate amplitude [C _1 r(— q, if)C]\ where C = 7274 and the 
index t denotes a matrix transpose. The pion decay constant /„. can be expressed as the loop 
integral 

UPp&ij = (0|<2-^759KjOP)} 

« l5ltl S(q + )r n (q;P)S(qj. (13) 
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For a chiral partner of the pion, a meson, we take the BS amplitude with only one covariant 

T a (q;P) = E a (P 2 )f (q 2 ). (14) 

At T = the mass-shell condition for a meson as a qq bound state of the BSE is equivalent 
to the appearance of a pole in the qq scattering amplitude as a function of P 2 . At T / 
in the Matsubara formalism, the 0(4) symmetry is broken by the heat bath and we have 
P — > (O m , P) where Q m = 2mixT . Bound states and the poles they generate in propagators 
may be investigated through polarization tensors, correlators or Bethe-Salpeter eigenvalues. 
This pole structure is characterized by information at discrete points Q m on the imaginary 
energy axis and at a continuum of 3-momenta. One may search for poles as a function of P 2 
thus identifying the so - called spatial or screening masses for each Matsubara mode. These 
serve as one particular characterization of the propagator and the T > bound states. 

In the present context, the eigenvalues of the BSE become A(P 2 ) — > A(f2^, P 2 ; T). The 
temporal meson masses identified by zeros of 1 — A(f2 2 ,0;T) will differ in general from the 
spatial masses identified by zeros of 1 — A(0, P 2 ; T). They are however identical at T = and 
an approximate degeneracy can be expected to extend over the finite T domain, where the 0(4) 
symmetry is not strongly broken. 

The general form of the finite T pion BS amplitude allowed by the separable model is 

r w (?„; P m ) = 75 (^(P*) + 74 Sl m F^P 2 m ) + 7 • PF^{P 2 m )) f (q 2 n ) . (15) 

The separable BSE becomes a 3 x 3 matrix eigenvalue problem with a kernel that is a gener- 
alization of Eq. pijl . In the limit fl m — > 0, as is required for the spatial mode of interest here, 
the amplitude F n = £l m Fn is trivially zero. 

The BS amplitude for the a at finite temperature has only one covariant. 

The impulse approximation for the <77T7r vertex, after the extension to T > for spatial 
modes characterized by Q = (0, Q) for the a and P = (0, P) for the relative tttt momentum, 
takes the following form: 

g<„rA T ) = -2A c T^tr s J ^^r 7r (A;„ + ;-P + )5'(g n+ _)r CT (g„ + ;Q) 

xS(q n++ )r n (k n _;P_) S(q n _). (16) 
The corresponding width of the a meson is equal to 



3 ^(T) 4M|(P) 
" l ] 2 16nM a V M2(T)' 1 } 



4 Numerical analysis and conclusions 

The results for tt, a and p mesons as well as related quantities are calculated with the parameters 
m = 6.3 MeV, b = 730 MeV, a = 0.52, A = 698 MeV, A 1 = 1.78 GeV, D A 2 = 223, 
Dih.\ = 137. This set allows one to reproduce the experimental data for = 140 MeV, 
f n = 93 MeV and M p = 770 MeV and obtain the reasonable result for the p — > tttt decay width 
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= 228 MeV (Vf?^ = 151 MeV). For the quark condensate we have (qq)° = (0.212 GeV) 3 
and the quark mass function at zero momentum equals m(0) = 484 MeV. 

We use this set of parameters to calculate physical observables of the scalar meson. As a 
result, we obtain M a = 762 MeV and Y a ^ = 456 MeV. 

In this work, we consider the T-dependence of the quark mass function m(p), f n and the 
spatial masses in the 7r and a channels. 

The T-dependence of m(0) , m (0) and f T is displayed in Fig. HJ ^o(O) corresponds to the 
solution of DSE in the chiral limit m = 0. The temperature DSE equation contains three 
functions A(p^,T), C(p^,T) and B(p^,T). The numerical calculations show that the solutions 
A(p 2 n , T) and C(p^, T) as functions of temperature practically coincide in the region < T < 200 
MeV. The temperature Td = 123 MeV(quark deconfmement temperature) presented in Fig. Q 
is determined by the properties of Green quark function in the DSE equation. It correspods to 
the appearance of the poles at real momentum in the lowest Matsubara quark mode. 

The critical temperature T c is that when the quark condensate disappears and chiral sym- 
metry is restored. 

In the chiral limit the pion is massless below T c and its mass increases above T c ; the a 
meson mass drops to zero at T c and above this temperature it is degenerated with M n . When 
m 7^ M n (T) is seen to be only weakly T-dependent up to near T c , where a sharp rise begins, 
as displayed in Fig. El These qualitative features of the response of the pion mode with T agree 
with the results deduced from the DSE in [T2]-|14j. 

The temperature dependence of the width r o ._ >W7r is also shown in Fig. El 
In conclusion, we have studied the ir and a meson properties and finite temperature in the 
framework of the nonlocal quark model which was used in jU] . The temperature dependence of 
these quantities shows that the restoration of chiral symmetry, M a ~ M n , occurs at the critical 
temperature T c = 130 MeV. 
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Figure 1: Temperature dependence of the quark mass function m(0) (solid line). Dashed line 
corresponds to the quark mass function in chiral limit, m (0). Doted line shows the temperature 
dependence of the decay constant f n . 



[5] T. Hatsuda and T. Kunihiro, Phys. Rept. 247, 221 (1994). 

[6] Yu. L. Kalinovsky, L. Miinchow and T. Towmasjan, Phys. Lett. B 283, 367 (1992); 

D. Ebert, Yu. L. Kalinovsky, L. Miinchow and M. K. Volkov, Intern. J. Mod. Phys. A 8, 
1295 (1993); 

D. Ebert, Yu. L. Kalinovsky and M. K. Volkov, Phys. Lett. B 301, 231 (1993). D. Blaschke, 
Yu. L. Kalinovsky, P. Petrow, S.M. Schmidt, M. Jaminon and B. Van den Bossche, Nucl. 
Phys. A592, 561 (1995). 

[7] P. Costa and M. C. Ruivo, Europhys. Lett. 60(3), 356 (2002). P. Costa, M. C. Ruivo 
and Yu. L. Kalinovsky, Phys. Lett. B 560, 171 (2003); Phys. Lett. B 577, 129 (2003); 
Phys.Rev. C70, 048202, (2004), P. Costa, M. C. Ruivo, C. A. de Sousa, Yu.L. Kalinovsky, 
Phys.Rev.D70, 116013, (2004); Phys.Rev. D71, 116002, (2005). 

[8] J.T. Lenaghan, D.H. Rischke and J. Schaffner-Bielich, Phys. Rev. D 62, 085008 (2000). 

[9] D. Blaschke, G. Burau, Yu. L. Kalinovsky, P. Maris, PC. Tandy, Int. J. Mod. Phys. A16, 
2267, (2001). 

[10] A. E. Radzhabov and M. K. Volkov, Eur. Phys. J. A 19, 139 (2004). 



6 



0,8 




0,00 0,02 0,04 0,06 0,08 0,10 0,12 0,14 0,16 0,18 

T [GeV] 



Figure 2: Temperature dependence of M n , M a and Y a . 
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